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ABSTRACT: 

This study of propellant fracture has been subdivided into three 
parts: (l) pertinent elastic fracture developments, (2) viscoelastic 

fracture, and (3) dewetting. 

In (l) a fracture criterion is developed for stress fields producing 
mixed mode crack propagation by drawing a comparison with fracture in 
uncracked geometries. Such mixed fields are prevalent for cracks in 
rocket motors with star shaped cavities. 

In (2) a means for predicting crack growth and velocity has been 
developed for linear viscoelastic meterials strictly on the basis of 
stress intensity factors. A fractiire characterization of critical stress 
intensity factor for a CTPB propellant has been completed, which facilitates 
the application of the theory. 

In (3) extensive dewetting e3q)eriments have been completed which 
have confirmed prior theoretical developments of a quasielastic dewetting 
theory. Not only has the concept of a stress criterion for dewetting 
been substantiated, but the form of the surface predicted by the theory 
has been verified. y 



FOREWORD 



This report summarizes the results of work on a research project 
entitled Propellant Fracture Studies, which is funded by the Naval 
Weapons Center, China Lake, California. 

The report is prepared in three sections. The first section 
contains ideas and observations concerning elastic fracture in combined 
symmetric and skew ssnnmetric stress fields. Validation of the theory 
is made with gldssy pol3rmers, and the work was performed as a precursor 
to formulating similar results for viscoelastic materials. 

Section 2 deals directly with viscoelastic crack propagation in 
S3nnmetric modes and develops a much more simplified theory than exists 
in the literature for handling viscoelastic fracture . The author would 
like to acknowledge many discussions of this material with LCDR C. 
Hertzler. 

Section 3 deals with dewetting, which is the phenomenon that makes 
propellant different from ordinary viscoelastic materials. Extensive 
studies of dewetting have been reported previously on this contract 
and this section continues the report of those studies. Acknowledgement 
is given to LCDR J. E Wood, who has contributed some of the results of 
his PhD thesis work, which is on this topic, for this report. 
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e , e e 
a ,h ,c = 


elastic boundary quantities 


V , V V 

a ,b ,c , = 


viscoelastic boundary quantities 
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strip width 




definition of boundary 


c = 


crack length 


E 


modvilus 


f,g>F»G,| = 


general functions 


G 

h 

> ^2 5 <^2 ~ 


shear modvilus 
heaviside step function 
stress invariants 


K 


bulk modulus 


^ = 


stress intensity factor 


r = 


polar coordinates 
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time of arrival of crack 


t = 


sheet thickness or time 


U. 
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displacement 


u = 


strain energy 


x,y 


cartesian coordinates 



V 



Greek Symbols 



a 

P 



e. . 

la 



\ 

V 

$ 




a, a ,a 
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SECTION 1 



SOME OBSERVATIONS ON ELASTIC FRACTURE UNDER 
COMBINED LOADING 



INTRODUCTION 

Mich of the traditional theory of fracture has been cast in terms 
of failure surfaces in principal stress space, which have been given 
physically descriptive names such as maximum principal stress, distort! onal 
energy, etc, and these theories have been applied to bodies that are 
assumed to contain no dominant flaws. On the other hand, in more recent 
times an equally strong theory in the form of "fracture mechanics" has 
been developed upon the premise that a crack already exists within the 
body. As a result of these two developments, fracture analysis is 
compartmentalized into problems of one type or the other. 

The two approaches have been further dissimiliar in the types of 
stress fields typically treated. While the use of combined stress 
states is very common in the instance of failure surface theories, it 
is not at ail. common to apply fracture mechanics to combined loading 
conditions, and in fact a general criterion is not available where more 
than one mode of fracture is present. In this r^ort an attempt is made 
to bridge between the two for elastic materials, linking them more 
closely together and providing a more unified perspective. By bringing 
both of these theories together, a general fracttare criterion has 
emerged for combined loading and mixed mode crack propagation. 
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Work on the problem of crack propagation in sheets under combined loads 

was undertaken by Erdogan and Sih^^^, who proposed an explanation of crack 

trajectories in plane two dimensional stress fields based upon experimental 

fol 

results obtained with FMM (polymethylmethacrylate). Later Panasyuk, et al^^^"^ 
Independently proposed similar explanations based upon their work with glass. 

[3] 

Cottrell , has examined the energy released from a system when a crack 
extends radially along an arbitrary trajectory. He found that the direction 
of propagation hypothesized by the above investigators corresponds to the 
direction of maximum energy release, which is consistent with the general 
fracture theory proposed by Griffith. 

In the experimental studies that were conducted relative to the before 
mentioned investigations, the stress conditions which precipitated crack 
instability were only of secondary interest; consequently observations are 
offered herein to complete the picture by describing the local conditions of 
fracture at a crack tip in an arbitrary two dimensional stress field by relating 
them to fracture conditions in unflawed geometries. This leads quite naturally 
to a fracture criterion representable in space. From this development 

further observations can be made regarding the changing nature of the stress 
field in the vicini.ty of a crack as it continues to propagate under combined 
loading. Verification of the theory has been demonstrated on PMM for two 
dimensional plane geometries. 

Interest in the combined loading fracture conditions arises for solid propellani 
rocket grains since many of the cracks in rocket motors axe subjected to both 
symmetric and anti symmetric loads. 
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COREEIATION OF FRACTURE CRITERIA 



For the unflawed case, i.e., for specimens without cracks, fracture 
in raultiaxial stress states is describable in terms of a surface in 
principal stress space. Over the years, it has been found by experience 
that the continuous siirfaces are quadric surfaces and are representable 
by second order equations in stress. A very general representation of 
such surfaces can be made in terms of the first two stress invariants 
I^ and Igj however it is physically more meaningful to use I^ and Jg, 
the second deviatoric stress invariant. The equation of the surface 
in functional form is 

f (I^^, Jg) = 0 (1-1) 

2 

Where l^ is selected for dimensional homogeniety. 

Mathematical expressions such as equation (l-l)' represent the manner 
in which a material reacts fracturewise to applied loads and theoretically will 
predict fracture regardless of the geometry or character of the stress 
state. A crack in a body introduces a geometrical change, but it 
should not influence the fundamental fracture behavior. Therefore a 
similar function should exist, which is related to this one, and which 
describes fracture of cracked bodies. 

[41 

For geometries containing cracks, Barenblatt has proposed that 
the local distribution of cohesive forces over the crack surface at 
fracture is a material property, independent of external loading 
and body geometry; thus for any body constructed of a given material 
there exists a function of stress intensity factors which describes 
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elastic fracture conditions. For two dimensionea stress states limited 
to the opening and sliding mode, his criterion reduces to 






( 1 . 2 ) 



where 

are due to the cohesive forces at the crack tip. The total stress 
intensity factor is obtained by adding the cohesive force contribution 
to and K^, which result from the external forces. 



and 



are the portion of the stress intensity factors which 



* >^a >=2t = * '^a • (1-3) 



By 



Barenblatt 




hypothesis of no singularity at the crack tip 



= -K^a = -'^a 

Thus from equations (1.2) and (l.U) 

§ (Kj_, K^) = 0 (1.5) 

To relate f in equation (l.l) to $ in equation (1«5), ^ 2 ^ will 

be expressed in terms of invariants. The stress field in the vicinity 
of a crack, when cohesive forces are not considered , may be written as 

'^io = h ^ij ^ ^ij 

From this equation, one can solve for and in terms of the two 
in-plane principle stresses or in terms of the two stress invariants, 
both of which can be used to relate fracture of cracked geometries to 
fracture of uncracked geometries. Continuing in terms of invariants, 
equation (1.^ can be vised to formulate and J 2 . 



V*' 



4 



(1.7a) 



r I. 



2 



(Kj^, Kg, 0) 



1 




2 (Kj^, Kg, e) 



(1.7b) 



Solving for K^ and Kg and inserting them into (1.5) produces a description 
of the local conditions governing fractvire at a pre-existant crack 



This describes the conditions of crack propagation in some direction, 

9^, of the maximum tensile stress, as referenced from the original crack 
orientation. It disregards the cohesive forces and expresses the 
fracture condition in terms of the externally applied loads only^ just 
as the fracture surface does. On the other hand, the success of the 
fracture surface concept over a wide variety of stress states and body 
shapes makes evident the fact that the fracture behavior of an elastic 
material is independent of the geometry of the body of interest, and 
it will not be altered in a fundamental way by the introduction of 
different geometric factors , such as cracks . It seems , therefore , 
appropriate to offer as a hypothesis that, when only external forces 
are considered, "Fracture will occur in cracked geometries at the same 
local stress conditions which produce fracture in uncracked geometries . " 
In other words, the crack will extend when the invariants of stress on 
the trajectory satisfy equation(U) and equation (1.8). This can be 
expressed mathematically as 




( 1 . 8 ) 



$ = f 



(1.9) 
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In other words, if the fracture surface f is known, the fracture 

condition, $, for bodies containing cracks may be obtained by replacing 
2 2 

Jg and with r Jg and rl^ . From this a general fracture criterion 
can be generated in terms of and Kg • 



EXAMPLE PLEXIGLASS 

Figure 1-1 shows the uncracked multiaxial fracture behavior of 

T6 1 

PMM obtained by Thorkildsen and Olszewski^ where it has been 
plotted in two dimensional principal stress space. A more precise 
reduction of this data in terms of invariants makes possible the 
correlation with results obtained by Erdogan and Sih^^^ on cracked 
specimens made from the same material. 

Stress invariants, I^ and Ig, are defined as: 



I^ = + CTg + 



= -(o^ Og + Og 03 + a ) 



(1.10a) 

(1.10b) 



A plot was made in Figure 1*2 of Ig vs I^ for the experimental data, 
which appears to yield a clearly defined linear relationship . The 
data was fitted by least squares to a linear equation, which resulted 
in 



= -2.93I2 + 0.960^^ 



( 1 . 11 ) 



where is the failure stress in uniaxial tension. This expression 
in terras of invariants is postulated to be the description of a 
failure surface for all possible tensile stress states. Instead 
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FIG. I“l FRACTURE STRESS FOR POLY Imethylmethocrylate) 
(Thorkildsen and Olszewski^®^ ) 
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FIG. 1-2 LEAST SQUARES FIT TO FRACTURE 
DATA FOR PMM 
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of leaving the final expression in terms of which has no particular 
physical interpretation, it is preferred to \ise J^, which is related 
to and by 




( 1 . 12 ) 



Ils-ing (1 -ii() and ( 1 . 12 ) an expression for the failure surface of PNIM is 
found to be of the form, 



f(ll^, Jg) = J2 + 0.0071^ - Const = 0 



(1.13) 



As can be seen, there is a very slight dependency on hydrostatic 
pressure, but shear is the predominant factor for fracture in tension 
of PMM. 

As stated previously, it is proposed that this criterion also 
governs the onset of crack propagation for bodies containing residual 
cracks in the same type of material. In other words, for bodies 
containing a crack, fracture will occur by an incremental extension 
of the crack along the path normal to the greatest tension when the 
stresses on that trajectory combine to satisfy the fracture criterion 
found for uncracked specimens . ^is requires that the direction of 
propagation, 0 ^, be known a priori; however this information can 
be calculated from the work of Erdogan and Sih. By equations (1.8) ^ (1.9) 
and (1.13), propagation of a crack is given by 

rJ2 + .007 rl^^ - Const =0 (l-l^) 
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To apply this criterion, the crack tip stress field is used. 

The coordinate description is standard, and the equations are classical 

[7] 

as given by Sih and Liebowitz . 



a = 
r 



^-T Cos ^ [(3 - Cos 9)K^ + (3 Cos 0 - 1)K2 Tan (l.l5a) 



2i2ry- 



^0 = 



2(2r) 



^ Cos I [(1 + Cos 0)K^ - 3K2 Sin 0] 

•r* 1 2 



(1.15b) 



\ 0 " 



2(2r) 



^ Cos I [K^Sin0 + K^{3 Cos 0 - l)] 



(1.15c) 



in polar coordinates is given by, 



jg = 5 - Ve * '^e * 3^6^ 



(1.16) 



which in terms of stress intensity factors from eqiiations (1.15) give 
r J2 = I [7 + 4 Cos0 - 3 Cos 20] + [3/2 Sin 20 - Sin 

+ I [7 - Cos 0 - I Sin^ 0]K^} (1.17) 

From equations (H5)> can also be readily computed, 

r = 2 [K^ Cos^ I - Sin0 + K^Sin^ |] (I.18) 

Combining equations (1.17) and (1 for a surface of the general form 
that was exhibited by PMM in equation (l.l^j we have 

r(J2 + al^) = I {[| (7+4 Cos0 - 3 Cos 20) + 6a (1 + Cos 0)] 

.+ [ (|- Sin 20 - Sin0) - 12 a Sin0] K^K2 (1.19) 

+ [| (7 - Cos0 - I Sin^0) + 6a(l - Cos0)] K^} 
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The constant in(l.3i+) is found hy examining the fracture condition 
of symmetric loading, for which = 0 and 9^ = 0. From (l-19) 



[r(J2 + = Const = |(l+12a) 



( 1 . 20 ) 



!Hius the criterion of fracture for arbitrary combined loads is given 
by 

(1 + 12a)K^^ = [(^ + 6a) + (| + 6a) Cos0 - | Cos 2e] 

+ [I Sin 29 - (l+12a) Sin9]lC|_K2 (1.21) 

+ [(| + 6a) - (| + 6a) Cos9 - ^ Sin^9] 
where for PMM, a = 0.007. 

To limit check the skew-symmetric fracture condition, = 0, 

Kg ■* Kg^ and Erdogan and Sih's theory and measurements give 9^ = -70.5°. 
Equation (l*2l) gives the critical Kg value in terms of the critical 
Kj^ value. 

Substantiation of this prediction from the theory is provided by the 
experimental results of Erdogan and Sih, who found for PMM, K^^^ = h'J2 
lb in and Kg^ = 422.1 lb in Using this value of K^^ in (1.22), 

and a = .007 from (l*l4)the theory predicts Kg^ = 421.5 lb in 
TOiis very close correlation provides meaningful substantiation to the 
theory. 
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Further evaluation of the theory can be obtained by plotting 
equation (1>2L) for a variety of combined loading conditions. For any 
given angle of propagation, 9^, the fracture equation (1-2L) is an 
inclined ellipse in K^, Kg space. However in K^, Kg space, each 
set of coordinates dictates a different loading condition, which 
requires the angle of propagation to be different, and the coefficients 
in the fracture equation vary, deviating from an ellipse. A plot of 
the equation is given in figure 1-3, where it is compared with the data 
of Erdogan and Sih for a variety of combined loading conditions 
obtained by using the geometry of a crack inclined to the direction 
of applied loads for several different crack orientation angles, p. 

As can be seen, the theory fits the data quite well for strongly 
skew-symmetric loadings, but deviates somewhat in some of the loadings 
that are predominantly symmetric. Of course for those loadings that 
are purely symmetric, the criterion matches identically due to the 
fact that the governing equation was matched to this condition. 



SUBSEQUENT BEHAVIOR 



After the crack has begun to run under combined loads, the 
geometry will appear as depicted in figure 1% for small extensions. 

A solution of this problem has been obtained by Anders s on ^ wherein 

he found that the stress intensity factors K_ and KU are 

/N r ^ ^ 



[4t^] ^ 



T 



(1.23) 
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FIG. 1-3 K, VERSUS Kg AT THE BEGINNING OF 
CRACK EXTENSION IN A CRACKED PLATE 
UNDER PLANE LOADING. 






FIG. 1-4 GEOMETRY OF CRACK AFTER EXTENSION 
UNDER COMBINED LOADS. 
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However before the crack extended, the stress intensity factors were 
those of an inclined crack in an infinite sheet. 

2 

= <J^ Sin 3 ^ ~ ^ ^ 

When these expressions are substituted into (1.23) it reduces to 

[K^ Cos 0Q - Kg Sin 0^] 

[-Kj^ Sin 0Q - Kg Cos 0^3 

This set of equations suggests that, in effect, propagation can be 
considered as a rotation of the K^Kg axis system, and from the 
experimental evidence available, the rotation is such that it reduces 
the value of Kg. In other words the crack, operating under the 
observed laws of behavior mentioned previously, will propagate in 
the direction of the maximum tension £ind this will reduce the value 
of Kg. Furthermore by extrapolation, subsequent increments of 
propagation will be eilong the direction of instantaneous maximum 
tension, which will continue to reduce Kg, driving it to zero so 
that ultimately propagation continues in the direction normal to 
the maximum tensile field stress of the body. 

To illustrate this, consider figure 1-4 as an example, where i = 0 
and p = 70°. The angle of maximum tension is 0^ = -30°, or -tt/6. 

The corresponding point on the fracture envelope in K^^, Kg space from 
(l.2l)and (1.24) is K^^ = 442 lb in and Kg = 155 lb in~^^^. After 
propagating an infinitesimal amovint, X, the new stress intensity 




(1.24) 



(1.25) 
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factors from (1.25) are - 44?. 4 lb and Kg = 84.3 18 

Thus- in the first infinitesimal increment of propagation, Kg is 
reduced by 46'jS. This change in Kg, which produces a change in the 
loading that the crack sees, causes a change in the stress field. 

These changes reduce Kg further, and this interaction continues until 
Kg is nonexistent and the crack is propagating in a symmetrically 
loaded condition. 

Another interesting way of looking at equation (1.25) is to take) 
the square root of the sum of squares. This results in 

(1.26) 

In Kj^, Kg space the square root of the sum of squares is the length 
of a positioh vector. Since ^(6^) ^ Ij this equation shows that the 
position vector to the fracture envelope in K space is becoming shorter, 
or that the "resviltant" stress intensity factor is growing smaller as 
the crack propagates. If it decreases by a sufficient amount, the 
crack shovild stop. The maximum reduction for the first increment 
occurs at the largest 0^, which is produced by pure shear loading. 

In that instance 0^ = -70. 5° » for which ^(6^) = O.85. This turns 
out to be enough to stop propagation of the crack after it has begun. 

In the experimental resvilts of Erdogan and Sih on PMM in pure shear, 
the crack subjected to pure shear stopped after propagating a short 
distance; however, this was the only crack geometry studied for which 
the crack stopped close to where it had begun to propagate. 
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These observations have been discussed within the context of 
experimental data on PMM and glass. Similar data on inclined cracks 
in plates are not available for metals. If they were, the validity 
of the observations could be checked further and the influence of 
plastic yielding could also be studied; however preparation of metallic 
sheet specimens with real cracks of arbitrary angles of inclination 
is quite difficult. Meanwhile the available evidence on two elastic 
materials appears to substantiate a correlation between the two 
failure domains . 



SUMt-lARY 

Observations have been made which link fracttire criteria in 
unflawed specimens to crack instability in bodies containing residual 
cracks, and from this a criterion for fracture \mder arbitrary loading 
is postulated in terms ofanenvd.ope in stress intensity factor space. 
Furthermore an explanation is offered for crack behavior subsequent 
to the onset of propagation which is based upon the concepts generated 
concerning the svirface in stress intensity factor space. Experimental 
results for two materials, glass and polymethylmethoerylate , verify 
the theory for two dimensional, combined loadings involving the 
opening mode and the hiding mode of crack propagation. 
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SECTION 2 



VISCOELASTIC FRACTURE 



SECTION 2 



VISCOELASTIC FRACTURE 

Introduction 

In the last three years viscoelastic fracture theory has been 
advanced along much more fruitful avenues than those which the previous 
efforts produced. By applying principles of conservation of energy at the 
crack tip, Knauss^^*^has shown that a differential equation governing 
crack length history can be generated. Several of his works continue this 
development, £ind in each case an unknown length parameter manifests itself 
in the differential equation. The parameter, which must be determined 
experimentally, seems to be inherent to the method of solution. 

Since this parameter is a bother to deal with and since it is 
artificial in the sense that there is no physical basis for its existence, 
it was felt that a worthwhile contribution could be made by devising a 
means of formulating viscoelastic fracture without the parameter. The 
means by which this was finally done was through the viscoelastic correspond- 
ence principle. 

The classical correspondence principle normally applied in the 
solution of boundary value problems in the linear, quasi-static theory of 
viscoelastic, does not apply for crack propagation problems. In that 
instance new boiindary is being created as the crack propagates, and the 
conditions on the boundary are changing in type; i.e. they are changing 
from displacement to stress bovindary conditions and vice versa. In this 
kind of problem there will be points of the boundary where only partial 
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histories of the boundary conditions are known; consequently the integral 

Laplace Transform cannot be computed, and the correspondence principle is 

not applicable. 

fQ] 

Graham has extended the correspondence principle to encompass 
problems of this type. (An explanation of his development is given in 
Appendix A.) By using Graham's method of solution, two very large 
categories of problems have been solved. Linking these solutions with 
Knauss' experimental findings and theoretical concepts has lead to a 
method of viscoelastic fracture analysis based upon stress intensity factors. 
The method parallels elastic fracture theory very closely and combines the 
work of other investigators in such a way that the length petrameter is 
no longer needed. 

Stress Analysis 

One very common class of crack problems for which the extended 
correspondence principle is applicable is shown in Figure 2-1, where 
the lower portion of the body has been removed and replaced with 

r 27 i 

appropriate boundary conditions. The two solutions by Graham in 
the literature that use the extended correspondence principle are for 
bodies of infinite extent. However the method is applicable to finite 
geometries. In the case of stress boxmdary conditions, like those 
pictured in Figure 2-1, the application is quite sinple. In the case of 
displacement boundary conditions, which is treated later, the bovindary 
conditions have to be recast into an acceptable form; therefore 
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Figure 2-1 Definition of boundaries for the 
Correspondence Principle 



the boundaries may be finite or infinite. They are loaded with 
normal stresses and are free of shear stress. The stress distributions 
are symmetric so that displacements on the imaginary extension in front 
of the crack are zero. The boundary conditions vary with time 
according to the following description: 



Q 

CQ 

n 

o 


on B 




= B(x,t) 


on Bj^(t) 


(2.1) 


u = 0 
n 


on 





These conditions correspond to column 2 of the matrix given in equation 
(A.I) of Appendix A. 

The extension of the correspondence principle by Graham permits the 
construction of the viscoelastic solution for a moving crack from an 
elastic solution for a moving crack. This in turn may be generated from 
the stresses and displacements of a stationary crack. For an incompressible 
elastic body, the field quantities in the neighborhood of the crack tip 
are given by 



a.. =-^f. .(e) 



(2.2a) 
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The crack is now imagined to propagate slowly. For most elastomers 
and all propellants, the velocity of propagation is quite low, which 
makes the inertial terms in the equations of motion small ajid neglectable. 
Therefore the stress and displacement fields are the same as those given 
in equation (2) except that the coordinate system to which the above 
quantities are referenced is now moving since it is attached to the 
crack tip. The dynamic, elastic solution is 



It is a simple matter to reference these quantities to an inertial 
coordinate system. As stated previously the loading and geometry is 
taken such that only the opening mode of propagation is present; there- 
fore 




(2.3a) 




(2.3b) 



r(t) cos0(t)+c(t) = X 



(2.U) 



y 




c(t) 



Figure 2-2 Inertial Coordinate System 



Equations (2.3) for the dynamic elastic problem became 






K^(t) 



F,,[e(t)] 



V2tl(x-c(t)) 

To use the correspondence principle, the transform of equation ( 
is needed. 

e* r '' 



IJ 



/y2rr[x-c(t )] 



F,,[0(t)] 



By equations (A. 4) of Appendix A, 



V* 

CT. . = 

ij 



K®(t) 



* 



/\/2n[x-c(t )] 



F,,Ce(t)] 

1 J 



E -E* 
v-*v* 



\ ' ^ V " 



2[x-c(t ) 



Gje(t)] 



E -E* 

V -♦ V*- 



For elastic stress intensity factors in which there are no material 
properties, equation (2.6a) is readily inverted and yields 

oT .(t) = a? . (t) 

ij 10 

From this result, by comparison of the terms, 

I^(t) HK®(t) 

where the stress intensity factor for viscoelastic materials is 
defined to be the same as for the elastic case. Indeed when the 



(2.4a) 

(2.4b) 

2.4) 

(2.5a) 

(2.5b) 

(2.6a) 

(2.6b) 

(2.7) 

(2.8) 
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boundary conditions are all in terms of stress, the two are identical, 

which indicates that the definition chosen is consistent. In other 

words, we are here defining the viscoelastic stress intensity factor as 

F. [e (t)] 

that factor which when miiltiplied by gives the stress 

V2n(x-c(t)) 

components in the vicinity of the crack tip. This def ini tion 
possesses the appropriate form as one proceeds to the limit of elastic 
materials. Also the stress intensity factor, contains information 
concerning geometry and boundary conditions, as it does in the case of 
elastic materials. 

For many geometries fitting the boundary conditions of equation 
( 2 . 1 ), the elastic stress intensity factor is 

= \Ov/?rc ( 2 . 9 ) 

where \ is a geometric modifiying factor that can be independent of time, 
but is often a function of c. If a is the applied stress at the boundary, 
then 

I^(t) = K^(t) = \[c(t)la(t)v^;^^ (2.10) 

Fracture Postulate 

Fracture studies in recent je ars have led to the conclusion that 
flaws in materials will grow when under load. This continues until the 
flaw becomes of critical size, at which time the specimen undergoes 
unstable crack propagation and fracture. This has been found to be 
true for metals, especially in fatigue, and for nonmetals alike. 

Knauss^^*^^ has published theoretical developments of these ideas 
for viscoelastic materials and has shown that the important parameter 
of interest for rate sensitive materials is crack velocity. In the 
elastic case, the threshold of rapid crack propagation is often the 
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focus of interest; however in viscoelastic materials the crack is 
found to propagate for almost any load level. The velocities are 
infinitesimally small for low loads and increase as the body is more 
highly stressed. The question to be answered then becomes, not at what 
load level does the crack propagate but what is the crack velocity for 
any given load. As defined, the viscoelastic stress intensity factor 
contains information relative to the boundary conditions; therefore 
■we desire the relationship between stress intensity factor and crack 
velocity. 

This brings us to the fractxire postulate: 



In words this states that for every stress intensity factor there is a 
corresponding unique crack velocity. 



have duplicated these tests and done additional work on other types 
of loads . ) For all temperature ranges tested, the material was elastic 
in the sense that it did not relax; however the crack behavior was 
viscoelastic and exhibited a velocity that varied with loading. Knauss 
postulates that in this type of material when the strain at the crack 
tip is the same in two different bodies, the same crack velocity will 
be produced. Of course in a body that is behaving elasticaJ.ly so far 
as stress and strain are concerned, the same strain level is equivalent 
to the same stress intensity factor. Thus in the postulate of fracture, 
we are following the previous ideas of Knauss. Furthermore, it 
states that if the stress intensity factor history can be calculated. 



If 
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the crack velocity equation can be determined from the experimental 
characterization of the material. Once this is known, the crack 
history can be obtained by quadrature. This will be demonstrated later. 

Fractiire Characterization 

Fractiire characterization hsis been done for only a few materials, 
and it has normally been accomplished with an infinite strip containing 
a semi-infinite crack. It could just as well have been done, however, 
on the sheet with the internal crack that was described earlier in this 
section. 

To apply the correspondence principle, we appeal to the extended 
version in Appendix A. The geometry is shown in Figure 2-3> where the 
boundary is given a normal displacement in the vertical direction on 




Figure 2-3 Displacement Boundary Conditions 



The actual boundary conditions are given below. 



ctY = 0 



u"^ = 0 
n 

u^= 0 



on B^(t) 
on B^ (t ) 
on Bg (t ) 
on Bg (t ) 



( 2 . 11 ) 



^s ®3 



on 
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These must he redefined in terms that will satisfy the conditions of 
the extended correspondence principle. This can he accomplished with 



the following set: 



= A"^(x,t) 



V 

aj= 0 

V 

Uc = 0 



^ = \(x,t) 
"I = Us(x,t) 



on Bj^(t) + B 2 (t) 
on B 2 ^(t) 
on B 2 (t) 
onBj 
on B^ 



The normal displacements are specifiable over the crack surface without 
material properties because the far field boundary conditions are on 
displacements. This now allows the construction of an equivalent 
elastic problem, whose boundary conditions are 



u^ = A®(x,t) on Bj^(t) + B 2 (t) 

a® = 0 on Bj^(t) + B 2 (t) 

uj = Uj^(x,t) on B^ 

u| = Ug(x,t) on B^ 

Since the transforms are known for all time. 



e* 

u 

n 


= A®*(x,p) 


on 


B^(t) + B2(t) 




= 0 


on 


B^(t) + B2(t) 


n 


= u/(x,p) 


on 




e* 

u 

s 


= Ug*(x,p) 


on 





represent the boundary conditions to the corresponding viscoelastic 
problem. From Graham the viscoelastic solution is obtained from 



(A.4). 
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The stress intensity factors for the strip are developed in Appendix B. 
Using the results of equation (B.20a) in equation (2.5) for v = ^, 



“ij -Vr ® 






-I * 



//2n[x-c(t)] _ 



rij(o) 



G. (tt) 6 
1 



9 = 0 



X < c , 6 = 1 
o’ 

X > c^, 6=h(t-t^) 



By the correspondence principle 



(2.12a) 



(2.12b) 



V* 

j. . 





e (t) 
o' ' 


j 3 ^REL P 


LA/2lT[x-c(t)] _ 



F. .(0) 
ij 



( 2 . 13 a) 



y i[x -c(tJJ 

1 2 o' *\J n 



G.(tt) 6 
1 



Inverting the transforms, the viscoelastic solution is obtained. 






(2.13a) 



ij 



e (t) 
o 



_/y2ir[x-c(T)] 



dT 



Pij(o) 



uY(t) G.(n) 6 

l' 2 '\J TT l' 



(2.li<a) 

(2.14b) 



It is evident from these equations that now the simpler expressions 
are the displacements, since the boundary conditions are of that type. 
The stresses are now more conqilex due to the material properties that 
are contained in the stress intensity factor. 

For a step strain history, which is most usually the one that is 
applied. 






(2.15) 
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The stresses and displacements reduce to 



aT,(t) aI^{ * TTa^ f dr)F„(0) 



(2.16a) 



10 



\ V2n(x-c^) 



o [2 tt&-c(t)]] 



3/2 



10 



u^(t) = (2.16b) 

The stress intensity factor, by comparison with equation (2.4a), and 
by the previously given definition would be 

K, (t) =V-TKx-c(t)]e < — + ttJ ^ dr } (2.1?) 

^ V 3 °L/2rrrx-c o r2irrx-cfT)lV/2 
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This is a very complicated expression , and one not very conducive to data 
reduction from laboratory results. However, a very good approximation 
can be made , which makes the data reduction quite simple. 

A distance out in front of the crack in the strip ,the stress field 
is that of a strip without a crack. If the specimen is long enoiigh 
compared to the crack, a large portion of the body will be subject to 
the stress given in equation (B.I7). With that approximation, the stress 
intensity factor is 

K^(t) (2.18) 

where a(t) is the stress acting on the strip computed on the basis of 
the net cross-sectional area, i.e. with the crack length removed. In 
that case, from equation (2.5a) 



e* A /3b 
a . . = A/TT” 
ij 



g(t) 



A/2rr[x-c(t)] 



F. .(0) 



(2.19) 
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Because there axe now no material properties in the expression 

V* e* , 
a . . = a . . and 
ij id 












A/2rr[x-c(t)] 

Therefore the viscoelastic stress intensity factor is 



( 2 . 20 ) 






( 2 . 21 ) 



Since the displacements are fixed at the houndary, the stresses in 
the strip are relaxing, and iC (t) will he decreasing in time as the 
crack propagates. This can he shown from the viscoelastic stress- 
strain relationships obtainable from the elastic solution of Appendix B. 
Equation (B.I7) for an incompressible material gives 



= 1^11 e 
V 3 o 



(2.22) 



By the regular correspondence principle, the stress in the viscoelastic 
strip is expressed in the classical form. 

dc (t) 



= 5 lo- 

For a step strain history, equation (2.15), 

= I % 

and the viscoelastic stress intensity factor of equation (2. 21) can 
be expressed in an alternate form 



(2.23) 



(2.24) 



Ki(t) e^ 



(2.25) 



Either equation (2.25) or (2. 21) can be used for data reduction. 
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Data Reduction 



For each different level of e there will he a different crack 

o 

velocity history produced. Unfortunately it is not constant in time 
because the stress in equation (2.21) is not constant. Parenthetically 
it should be stated here that if a means for experimentally applying a 
constant stress boundary condition could be employed, then the resiilts 
of the first case of Appendix B could be used and nothing would vary 
with time in the strip. For a given level of stress the stress intensity 
factor would remain constant, since it is independent of crack length; 
consequently the crack velocity woiild be constant, and it would be a 
very desirable test. In fact it would be just like the elastic case 
for which this test has been used by Knauss for a number of years 

Returning to data reduction, the stress and crack length in the 
strip are recorded as a function of time, as shown in Figure 2-4. 





Figure 2-4 Representative Data From Strip Test 

From the crack history, the crack velocity history can be constructed by 
taking the slopes graphically (Fig. 2-5). A cross-plot can then be 
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made of stress, or equivalently stress intensity factor, versus crack 
velocity. This is the master curve, or fracture characterization 
cijrve shewn in Figure 2-5. Normally these tests will be conducted over 
a range of temperatures in order to obtain the whole curve. 





Figure 2-5 

Master Fracture Characterization Cvirve 

ri3i 

Hertzler has obtained such characterization data over a range 
of temperatures from 15°P to 85 °F for a CTEB propellant. He 
found that the data shifted very well to make a master curve, and the 
shift was exactly the same as that obtained by other investigators for 
the relaxation nodulus. His results are shown in Figure 2-6. Part 
of the scatter in the data is attributed to the fact that the approximate 
expression was used to reduce it. 

To our knowledge this is the first viscoelastic fracttire characteriza- 
tion of propellant. Some scatter in the data is caused by local inhomo- 
geneities in the material, which in turn cause the crack velocity to 
fluctuate as the tip proceeds through, or around, these regions. 
Nevertheless the master curve is well defined and exhibits no more 
scatter than is normally noted in propellant data of other varieties. 



30 




lO 

d 



o 

d 



tn 

d 

I 



q 

I 



m 



Q 

CVJ 

I 



iq 

CVJ 

I 



0 

ro 

1 



lO 

rd 



q 

I 



I 



31 



log CQx (in/min) 



even modulus data, which should evince the smallest amount of scatter 
of all. 

The master curve indicates two limit values that are important. 

The low velocity region gives evidence of a lower hound of stress 
intensity fantor for which there will be no crack growth. This figure 
of = 10 psi ^/in is a useful guide in making limit checks on behavior 
of cracks in rocket motors from -which their structural stability can be 
determined. Any loading below this limit should be safe, at least for 
the present time; however, the tendency of this characterization to change 
with time, as in a storage condition, is not known. Future testing 
could be directed toward this end; however it would require considerable 
amounts of propellant, since many rather large specimens are required. 

Specimens selected for this characterization were marginally small 
in size and were 6" x 1" x 1/4" and bonded to wooden tabs along the 
six inch dimension. Five different specimens were ru» for each of 
the six temperatures. Of course it would be best if the entire set of 
tests were duplicated at each time interval of the period selected. 

From this information a good description of fracture behavior over the 
desired shelf life could be obtained. 

In lieu of that data, one might assume on the basis of the results 
in Section 1 that the stress intensity factor -vri-ll degrade in the same 
manner as the ultimate strength. This quantity is routinely measured in 
biaxial endurance data and forms the basis for extrapolating to the 
storage life reqiiired. With this information, the behavior of cracks 
existing in rocket motors that are cvirrently stored could be predicted, 
and the reliability of the missile from this standpoint covild be assessed. 
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The other end of the spectrum is also of interest. It appears 
that if the data were extended in the high velocity range, the ultimate 
limit of crack velocity would be approximately 10 inches per minute. This 
is indeed slow and substantiates the earlier assumptions that neglected 
the inertial terms in the equilibrium equations. An ultimate limit 
strongly influences the crack history. For instance consider an infinite 
sheet with a central crack subjected to a step stress at infinity. 

From equation (2.IO) 



begin to grow. As it grows Kj^(t) increases, and in turn the crack 
propagates faster. This process continues until the crack reaches a 
limit velocity, at which time it will continue to propagate across the 
infinite sheet at constant speed, even though the stress intensity factor 
continues to increase. 

The time to reach the constant velocity portion of the history is 
a function of the magnitude of a^. If were large enough, the crack 
could begin its propagation in the constant velocity region, in which 
case the behavior would be elastic, much like most observed behavior with 
metals or glassy polymers. On the other hand, infinite time could be 
required if K(0) is below 10 psi >\yTn . However any X(0) above 10 psi \/in 
eventually will produce the constant -maximum velocity propagation for 
this kind of geometry, where K^(t) is continually increasing. 




(2.26) 



From Figure 2-6, if Kj^(o) is greater than 10 psi >s/inj the crack will 
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Solution Method 



The method for predicting the crack behavior in viscoelastic 
materials essentially follows the method proposed by Knauss^^^ for 
Solithane^with some variations for the relaxation behavior. For 
discussion purposes, consider a rocket motor with a bore crack, where 
it is desired to predict what the crack will do under a given loading 
program. 

If an analytic expression for the stress intensity factor of the 
elastic motor geometry can be obtained, then by the correspondence 
principle previously discussed, the viscoelastic stress intensity factor 
is obtednable. Suppose it has the form 

K^Ct) = l[c(t)la(t)vW (2 

From the fracture characterization curve similar to Figure 2-6, 

K^(t) = f(6) (2 

Inserting equation (2.28) into (2.27), a first order differential 
equation in crack length is produced, which when integrated will give 
the desired history. Comparison of that history with storage life, 
or burning rates under operational conditions, will determine the 
detrimental effect, if any, that the crack may contribute. 

Often times motor geometries are too complex to allow analytic 

expressions for the stress intensity factor. In this case numerical 

methods have been developed to find elastic stress intensity factors. 

r 19] 

Deverall and Lindsey have suggested some simple methods by which 
these can be obtained. Subsequently numerical transforms and inversions 
would be necessary to obtain the final solutions. Because this is an 



involved solution process, an approximate method is being developed and 
studied, which looks very promising. It is based upon an observation 
by Dr. Russell Westman of UCLA that the local stress in bodies without 
cracks can often be related to the stress intensity factor for the same 
body with a crack. Progress on this development will be given in a 
later report. 

Conclusions 

It is felt that the viscoelastic crack behavior for materials in 
which the crack propagates slowly is predictable from the concepts 
developed in this section. Limitations on loading and geometry, 
however, must be such that symmetric crack behavior is produced. Further 
investigations will be made to ascertain if the concepts developed in 
Section 1 for mixed mode propagation can be applied to viscoelastic 
fracture. Further evaluations of the validity of these concepts will 
be made using existing data from the literature and in generating new 
experimental crack propagations results where needed. 
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DEWETTING 



SECTION 3 



DEWETTING 



Introduction 

Propellant differs from classical viscoelastic materials through 

dewetting. Since dewetting has a pronounced influence on the stress- strain 

law for strain levels experienced in many rocket motors, it is necessary 

to investigate the influence of dewetting on fracture. 

r 173 

In previous reports Lindsey ' and co-workers have developed 

an isotropic theory for dewettable materials. From this development a 

means has been formulated by which stress analyses can be carried forth 

while accounting for dewetting. In that theoretical treatment it has been 

postulated that a stress criterion for dewettihg exists and that is is 

expressible in terms of a surface in principle stress space. Furthermore 

[l8] 

by application of Drucker's Fundamental Theorem for stable materials, 

[19] 

Wood has obtained the general expression that describes the dewetting 

criterion for arbitrary stress states 
2 

Jg + or I^ = Constant (3-l) 



In order to evaluate these results and the associated assumptions, 
experimental tests were developed and conducted for dewetting in uniaxial 



where 




2 
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and triaxial stress fields. 



The experiments employed a liquid dilatometer, both pressurized and 
unpressurized, and the pokerchip test . A brief synopsis of the experi- 
ments is presented here. At a later date the detailed treatment of the 
theory and experiments will be presented in LCDR J. E. Wood's PhD thesis. 
Dilatometer Tests 

A liquid- filled dilatometer using Dow-Corning DC-200 silicone fluid 
was employed for the measurement of volume change in uniaxial tensile 
test specimens. Load was sensed from strain gauges mounted on a flexure 
in line with the test specimen within the dilatometer chamber. Volume 
change was determined by visually measuring the fluid height in a capillary 
tube which permitted a minimum resolution of dilatometric strain of 0.01 per 
cent. In order to minimize the effect of thermally- Induced variations, 
fluid height readings were compared with readings from a reference chamber 
which duplicated the test chamber in all thermal characteristics. Maximum 
test duration was limited to thirty minutes. 

A number of different formulations of highly filled composite solid 
propellant were tested. However the most conplete set of tests was 
carried out on a highly filled PBAN propellant supplied by United Technology 
Corporation. All of the following data were obtained on that propellant. 
Constant Load Dilatometer Tests 

Constant load uniaxial tension tests were conducted in the dilatometer 
to investigate the time -dependent nature of the onset of dewetting under 
small constant loads. Figure 3-1 shows dilatation vs. time for a series 
of five tests of the UTC PBAN. The rapid increase in dilatation observed to 
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Figure 3-1 

DILATATION vs. TIME AT CONSTANT 
LOAD, UTC PBAN 




Figure 3-2 

POISSON'S RATIO vs. TIME AT CONSTANT 
LOAD, UTC PBAN 
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occur soon after the application of load for the three higher stresses 
gives a clear indication of dewetting. Because of the very sTnqn 
amounts of volume change produced in these tests at the lovjer stress 
levels, it is difficult to determine the presence or absence of dewetting 
from the figure. Further amplifying information is obtained from the 
variation in Poisson's ratio with time for the same specimens, shown in 
Fig. 3-2. In this plot the rapid decrease in u associated with dewetting 
is seen to occur immediately at 55.^ psi and 31.5 psi. The response for 
the three lower stresses shows that the drop in v is delayed in time and 
reduced in magnitude with decreasing stress. 

Qualitatively similar behavior was observed for all composite pro- 
pellant tested, including PBAN, PBAA, CTPB, and HTPB formulations. 

From these tests, it was concluded that although the onset of 
dewetting in the propellants tested does demonstrate some time dependence 
for sufficiently small values of stress, the small amoimt of dilatation which 
results maies evaluation of this dependence impractical in the dilatometer. 
However, a reasonably narrow zone of stress values exists above which dewetting 
effects are negligible, at least for loads of short duration. It is in this 
sense that a stress criterion for dewetting will be interpreted. 

Constant Strain Rate Dilatometer Tests . 

These tests were carried out to identify the stress criterion for the 
onset of dewetting in uniaxial tension and uniaxial tension with superimposed 
hydrostatic pressure at a strain rate of one per cent per minute, for which 
the material response prior to dewetting was quasi- elastic. 

Figure 3-3 shows the stress- strain-dilatation response for four dupli- 
cate tests of PBAN specimens. The dilatation increment plotted has the 
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Run Na 




Figure 3 — 3 

STRESS AND DILATATION vs. STRAIN-UTC PBAN CONSTANT STRAIN 

RATE TESTS 



initial linear response subtracted out, thus the intercept of the dilata- 
tion curve with the axis indicates the onset of dewetting. 

Following the above tests, pressurized dilatation tests were under- 
taken in 6in attempt to determine the onset of dewetting under a combined 
state of stress. Tests were run with superimposed pressures of 10, 20, 
and 40 psig. Maximum test pressure was limited to 40 psig because of 
excessive creep of the plexiglas walls of the dilataneter chamber experienced 
at higher pressures. More data scatter was encountered with these pressurized 
tests; however, no significant differences from the unpressurized tests 
either in the stress at onset of dewetting or in the character of the 
response could be detected. 

Pokerchip Triaxial Tension Test 

In order to investigate dewetting behavior in a stress state approxi- 
mating hydrostatic tension, a pokerchip apparatus was employed. For 
materials with Poisson's ratio near 0.5, the pokerchip triaxial tension 
test provides a state of approximately hydrostatic tension within the 
central region of the test specimen. Using the dewetting criterion previously 
discussed, it is foTond from the analytical solution for stresses within the 
pokerchip that the dewetting criterion is first satisfied at the center of 
the specimen, within the region of hydrostatic tension. Thus the poker- 
chip test provides a means of initiating dewetting under a known triaxial 
tension state. 

In order to interpret the results of the pokerchip test properly, 
it was necessary to determine the influence of dewetting on experimentally 
measurable quantities. To Investigate this question a computerized finite 
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element analysis of the pokerchip was performed. The Rohm and Haas 
AMG 032A structural analysis program for axisymmetric structures and 
loads was adapted to the IBM 36O-67 con5)uter at the Naval Postgraduate 
School. This program is designed to handle materials with Poisson's 
ratio at or near 0.5, and allows variation of material properties from 
element to element, thus it was suitable for the task. 

To test the program and establish convergence, the program was 
used to analyze a pokerchip specimen loaded with a uniform axial 
elongation. Results were compared with the approximate analytical 

C 2o3 

solution given by the LSWZ^ analysis. The geometry chosen was that 
which would be used in later testing - a disk 4" in diameter and 0.25" 
thick. Various modulus values representative of composite solid pro- 
pellant were investigated, with Poisson's ratios in the range 0.49 "fco 
0.499- The program was foxind to perform satisfactorily, giving close 
correlation with the LSNZ solution. The program was designed to have the 
capability of handling incoit^jressible material (u = .5), however, it was 
found that for the very stringent displacement boundary conditions imposed 
by the pokerchip test, convergence became less easy to achieve as u 
approached 0.5. Py utilizing the maximum allowable number of elements in 
the radial direction, it was possible to obtain good solutions for u = 0.499* 
Table 3-1 gives a comparison of stresses calculated by the program and by 
the LSWZ solution for this value. The comparison is quite good, diverging 
slightly only near the outer radius where the LSVE solution posseses 
singxilarities . 



42 



To model the effect of dewetting in the pokerchip test, it was 
considered that the initial effect of dewetting would be the alteration 
of material properties in a small region at the center of the specimen. 
For the finite element representation, a two-phase specimen was employed, 
having an inner core l/2" in diameter to which reduced values of E and u 
were assigned, while the properties in the rest of the pokerchip were 
unmodified. Figure 3-4 shows a cross-section of the finite element grid 
enployed. 

Results 

It was found that most aspects of the solution for the two-phase 

specimen showed only slight changes as properties within the core were 

altered. But a significant change was noted in the state of stress 

around the center of the specimen, especially within the core itself, 

for slight variations in Poisson's ratio in the core materieil. Figure 

3-5 shows the computed axial stress at the midplane versus radius for 

veirious values of u at the central core. For a reduction of u of 

core 

0.8^ (from .499 to .495) the stress at the center drops by approximately 
10^, with progressively greater reductions occurring as is further 

decreased. The indicated percentage change in the stress state at the 
center is at least an order of magnitude greater than that occurring in 
any of the customarily measured quantities such as total load or radial 
displacement on the periphery. 

On the basis of this analysis it appears that a pressure transducer 
located at the center of one of the platens of the test rig could 
easily detect a deviation in the local stress state due to a small 
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Figure 3-4 

FINITE ELEMENT REPRESENTATION OF POKERCHIP WITH 

DEWETTING 
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FINITE ELEMENT SOLUTION 
Poktrchip- radius = 2.0” 

thickness =0.25" 

E = 300 psi 
V = 0.499 psi 

radius = 0.25" 

E^ = 300 psi 

^ as shown 
Loading ^ axial displ. 0.00125" 
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Figure 3-5 

AXIAL STRESS vs. RADIUS In TWO-PHASE 

POKERCHIP 
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TABLE 3-1 Midplane Stresses in Pokerchip Test Specimen 



Specimen Diameter: 4.0", thickness 0.25" 

Material Properties; E = 300 psi, u = 0.499 
Loading: Uniform axial elongation - Az/ t = 0.005 



Radius 


Finite Element Sol. 


LSWZ Sol. 


Finite Element Sol. 


IiSWZ Sol • 


Inches 


a 

r 

PSI 


PSI 


PSI 


a 

z 

PSI 


a 

z 

PSI 


.0625 


74.5 


74.5 


74.1 


76.2 


75.4 


.1875 


73.8 


73.8 


73.6 


75.7 


75.0 


.3125 


72.8 


72.8 


72.5 


74.7 


73.9 


.4375 


71.1 


71.2 


70.9 


73.0 


72.3 


.5625 


69.1 


69.1 


68.8 


71.0 


70.1 


.6875 


66.3 


66.3 


66.0 


68.2 


67.5 


.8125 


63.1 


63.1 


62.9 


65.0 


64.2 


.9375 


59.1 


59.1 


59.0 


61.0 


60.4 


1.0625 


54.8 


54.8 


54.6 


56.7 


56.0 


1.1875 


49.5 


49.6 


49.6 


51.5 


51.0 


1.3125 


43.9 


44.0 


43.9 


45.9 


45.3 


1.4375 


37.4 


37.5 


37.6 


39.5 


39.0 


1.5625 


30.3 


30.4 


, 30.5 


32.4 


32.0 


1.6875 


22.5 


22.6 


22.8 


24.5 


24.2 


1.8125 


13.7 


14.2 


i 4.3 


16 .6 


15.8 


1.9375 


3.8 


3.5 


5.0 


5.1 


6.5 
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alteration of properties. While at present the extent to which 
dewetting in hydrostatic tension does alter the Poisson's ratio 
has not been determined, such an installation would appear to be at 
least as sensitive to Poisson's ratio as the dilatometer is for uni- 
axial testing, and should provide an equally good Indication of dewetting 
initiation. 

Another finite element problem was worked, which provided for a 
section of flexible boundary, just above the dewetted core, -where the 
pressure transducer would be mounted in the otherwise rigid grip. This 
analysis was conducted to determine if introduction of the pressure 
transducer affected the stress field in the specimen. Even by allowing 
the displacement at the boundary to be 200 times the expected displacement 
of the transducer diaphragm, the change in the stress field was imperceptable . 

Accordingly a pressxire transducer with a l/2" diameter contact area 

was imbedded in the center of one of the 4" diameter pokerchip platens, 

[21] 

similar in construction to that reported by Harbert . Figure 3-6 
shows the stress- strain plot obtained from one of the pokerchip tests run 
at an axial strain rate of 0.007/min. The dashed lines indicate the elastic 
solution for the specimen -with E = 520 psi and v = .495* At an axial strain 
of 0.0075 > the stress sensed by the pressure transducer at the center of the 
specimen shows a sharp reduction from the initial linear slope. The stress 
level at this point, 94 psi in this case, is taken as the point of de-wetting 
initiation. 

Pokerchip tests instrtimented with a central pressure transducer have 
been reported by Harbert. For his tests performed at a constant rate of 
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Elastic Solution 
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.02 .03 .04 

(overage) 



axial elongation, a sharp break in the pressure transducer reading was 
observed at a characteristic level of stress which was identified as 



the point of initial failure. From the analysis reported above, it 
appears that such a result may indicate only a slight reduction in 
load- carrying capability which may be at a stress well below the actual 
failure stress. 

Comaarison of Experimental Data with Theoretical Criterion 

The theoretically derived dewetting criterion referred to previously 
is 

2 

Jg + = Constant (3.1) 



is a material parameter dependent upon the change in the shear and 
bulk moduli of the material as dewetting initiates. From data for the 
unpressurized uniaxial tensile test runs shown in Fig. 3-3 > Q? was calcu- 
lated to be 0.0026, making the dewetting criterion for this material at 
the temperature and rate tested: 



The stresses measured at the point of dewetting in the various tests 
described are plotted in Fig. 3-7. The solid curve shows the theoretical 
dewetting criterion. 




Qt 




dG 



+ 0.0026 = 180 

2 

For the uniaxial tensile tests Jo = /3 and I, = cr - 3P 

2 u' 1 u 

For the pokerchip triaxial tests Jg = 0 and = 3^^^ . 
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UNIAXIAL TENSILE STRESS vs. DILATATIONAL STRESS at DEWETTING 

INITIATION, UTC PBAN 



While the limited amount and range of the data presented cannot 
alone confirm the validity of the stress criterion, it is satisfying 
to note the good correlation with the triaxial tension results. 
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SECTION 4 



CONCLUSION 



SECTION 4 



CONCLUSION 

Several aspects of polymer fractvire have been discussed in this 
report. A theory for crack instability under mixed mode conditions has 
been proposed, and to a limited extent verified, for a glassy polymer. 
Normally advanced concepts such as this are first attempted on the 
simpler elastic materials where there exists a higher probability of 
success. Then once a foothold has been established, it can be generalized 
and made applicable to viscoelastic materials. The developments and 
results in Section 1 constitute work of that type and further work will 
be directed toward generalizing it to be applicable to viscoelastic 
materials and propellants. 

Viscoelastic materials, being more conplex, do not have so advanced 
a fracture theory as do elastic materials! The evolution of theory lags 
that of elastic materials but is following along in a very similar pattern 
of development. The concept of stress intensity factor for viscoelastic 
materials is used very extensively in this development, much more so than 
in the current literature on the subject. The methods proposed by which 
materials may be characterized for fractvire and later analyzed for 
fracture are still hypothetical in nat\ire and must be evaluated 
experimentally. This work is presently being pursued. 

An alternate method of development based upon an energy approach as 
performed by Knauss and his coworkers, leads by a very complex analysis 
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to a dependence in the solution on an undetermined length parameter, 
which must he evaluated experimentally. The approach outlined in 
Section 2 avoids this problem by going directly to the stress intensity 
factor, and it becomes much more simple mathematically. 

Finally propellant fracture to be complete must incorporate dewetting 
effects. A theory to do this in quasielastic situations has emerged and 
has been experimentally evaluated for a limited number of stress states. 
Considerable work has been expended in developing this theory and LCDR 
J. E. Wood has completed a PhD research project in this area. His 
thesis is now in preparation and will treat this topic in detail. Copies 
will be forwarded when available. A brief synopsis of just a few of the 
IB suits are included in Section 3 to provide an indication of the type of 
work that has been done. 

An evaluation of what analyses can be performed at this time leads 
to the following conclusions: 

1. Predictions of crack instability for purely elastic materials 
can be made for two dimensional, through the thickness cracks 
under arbitrary loading. 

2. Predictions of crack history can be made for purely viscoelastic 
materials with no dewetting for two dimensional, symmetric 
geometries and loads. 

3. Dewetting can be incorporated into the stress analysis for 
quasielastic loading conditions such as cool-down. Study of this 
problem is presently underway and the associate problem of the 
temperature dependence of the dewetting criterion is also being 
studied. Measurements of dilatation over a range of temperatures 
has been completed. The step to incorporate dewetting into 
viscoelastic materials has not yet been made. 
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APPENDIX A 



EXTENDED CORPESPONDENCE PRINCIPLE 

Tq] 

The extension of the correspondence principle by Graham ^ to 
cover crack propagation problems is a key step in the development of 
viscoelastic fracture as it is done in this report. As a result, a 
synopsis of the central theme of the paper is included here; however 
the development does not parallel that of the paper but enlarges upon 
it. It is hoped that by reading both this appendix and the paper, a 
clear view of the basis for Section 2 of this report will be obtained. 
Definition of the problem 

The correspondence principle is clearly applicable whenever the 
tjrpe of boundary condition prescribed is the same at all points of the 
boundary. For mixed boundary value problems (i.e., problems for which 
different field quantities are prescribed over separate parts of the 
boundary) the method is still applicable provided the regions over 
which different types of boundary conditions are given do not vary 
with time. (We are, of course, assuming that the region occupied by the 
body does not vary with time . ) There remain those viscoelastic mixed 
boundary value problems where the regions, over which different types 
of boundary conditions are given, do vary with time. Particular 
examples are indentation and crack propagation problems. For problems of 
this type there will be points of the boundary at which only partial 
histories of some field quantities will be prescribed. When this is the 
case the transforms of these quantities are not directly obtainable and 
the classical correspondence principle is not applicable. 
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We desire the solution to the viscoelastic problem shown 



diagramatically in Figure A-1. 




FIGUEE A-1 

Viscoelastic Crack Problem 



In figtire A-1, the material below the line of crack propagation 
has been removed and subsequently replaced with the appropriate 
boundary conditions. a^(x,t) , b^(x,t) andc^(3qt) represent the elements 
of any column of the mactrix in A.l. 
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u^ and Ug denote normal and tangential components of the boundary 
displacements, and cr^ and represent simileir quantities for stress. 
Ilie boundary conditions for the problem are: 



a^(x, t) = A^ (x, t) 


on B 




b'^(x,t) = B'^(x,t) 


on B^(t) 


(A.2) 


c^(x, t) = 0 


on B2(t) 
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where B^(t) + BgCt) = B, B^(t) is monotonically increasing, and B 
is the boundary of a fixed region R. 

The solution is obtained by first solving an elastic problem 
with changing boundaries as shown in Figure A-2. 






e e 

ei^(x,t), u^(x,t) 



are solutions in R. 




FIGURE A-2 

Equivalent Elastic Problem 

0 

With the solution to this problem, the boundary values of b (x,t) 
on B2(t) are computed, and the values of c^(x,t) are computed for B^(t). 
With these values, now define 

b®(x,t) = B®(x,t) on 

c®(x,t) = k(E,u) C®(x,t) on 

The constraints on the problem must be such that the new boundary 
conditions B^(x,t) are independent of material properties. This is 
necessary for the classical correspondence principle to work in the next 
step of derivation. 

We can now define the same elastic problem in terms of new 
boundary conditions as shown in Figure A-3 



(A. 3) 
B 
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FIGURE A- 3 

Elastic Problem with New Boundary Conditions 
Because of the nature of the new boundary conditions, transforms can be 
computed for them since they are now defined for all time. That 
condition is shown in Figure A- 4. 



* * * 

CT. . , e. u. are solutions in R 




FIGURE A- 4 

Corresponding Elastic Problem 

This is the corresponding solution for the viscoelastic problem 
shown in Figure A- 5. 
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aT.(x,t), eT.(x,t), uy(x,t) are viscoelastic solutions in R 
1 j 1 j 1 



I 




FIGURE A- 5 

Viscoelastic Problem Equivalent to Figure A-1 



Therefore the correspondence principle is given by the 
classical relationships . 




(A. 4) 



This extended version of the correspondence principle is used 
to generate theoretical developments pertaining to viscoelastic fracture 
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Extension of Graham* s Paper 



An extension of Graham's development can he made by introducing 
a third boundary, B^j which is constant in time; so that now the 
unchanging body region, R, is bounded by 
B = Bj^(t) + B^it) + B^ 

Since B is constant, 

B^(t) + BgCt) = constant, 
as in Graham's paper. 

B^ can consist of an arbitrary number of segments, upon each of which 
the boundary conditions may be different; however upon none of 
the segments can the type of bovindary condition change with time. 

Such a difinition of the boundary conveniently facilitates the 
application of the correspondence principle to finite geometries. 

For instance consider the infinite strip geometry with displacement 
boundary conditions as shown in Figure A-6. 




Figure A-6 Infinite Strip Subjected to Displacement Boundeiry Conditions 
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The houndary conditions on are known for all time, and therefore 
their transforms can he computed. Now applying the original methods 
of Graham, the entire development goes through as before; however now 
a more versatile geometry is available for which the correspondence 
principle is applicable. The prime limitation upon generating solutions 
is the requirement that B®(x,t) in equation (A. 3) he independent of 
material properties. 
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APPENDIX B 



STRESS ANALYSIS OF AN INFINITE STRIP 
Viscoeleistic crack propagation has Been extensively studied by 

r 10 11 22 23 2UI 

Knauss and coworkers.^ * * > ->> j Much of the work has been 

done on an infinite strip with a semi infinite crack. A simplified, 
but oft times adequate, solution method is given in this appendix for 
the stress intensity factor of the strip. The elasticity theory 
required is very elementary, and it is presented here for con5)leteness 
in the development of the viscoelastic stress intensity factor. 

Elastic Stress Analysis of Strip 

(T 
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Figure B-1 Stress Boundary Conditions 
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T =0 

xy 



The Stress Boundary Conditions are specified 6is 

y = ^ 

No X = constant boundaries (infinite) 

Zero stress in the x direction at infinity. 

Consider a stress function of the form 

Obviously this function satisfies the biharmonic equation. The 
components of stress produced by this function are: 

= 0 



(B.l) 



a =0 

X 



dy = a 



xy 



(B.2) 



This stress field satisfies the boundary conditions and therefore 
constitutes the solution to the problem. The corresponding strain field 
is : 

e = - ^ a e=| c=-^a (B.3) 

If a semi-infinite crack is inserted in the strip, as shown in 
Figure B-2, the solution obtained above will be valid out in front of 

4 

the crack at some distance where the disturbance caused by the crack 
has decayed to zero. 



. i f I I I I I I I I I I I I 



^ ‘ I II M I 1 M I I 1 



Figure B-2 Infinite Strip with a Seminfinite Crack 
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To the fax left, the sheet is conipletely unloaded and stress free. To 
the far right it is under constant stress as if there were no crack, 
and in the center there is a transition region of unknown character 
between the two constant fields. If the crack advances by an amount Ac, 
the transition region moves ahead Ac. The constant stress field is 
reduced in dimension by Ac, and the unstressed region is increased by 
the same amount. The strain energy lost from the constant stress field 
by the crack moving ahead Ac is 

AU = i a e 2btAc (B.4) 

where t is the thickness of the sheet. In the limit of infinitesimal 
increments in crack length. 



= a e bt = ^ bt 
dc y y E 



(B.5) 



From Irwinfe development, which related stress intensity factor and 
strain energy release rate for cracks in plane stress, 



K, 



1 ^ 
t dc 



(B.6) 



Thus the stress intensity factor for the infinite strip is 

= oVb (B.7) 

It is of interest to note that this geometry is independent of 

crack length and therefore possesses a constant stress intensity factor 

as the crack propagates across the specimen. For this reason this 

geometry lends itself well to laboratory measurements. Mueller^ 

has obtained a solution of a pressurized crack in a semi- infinite strip, 

and from his results, one can determine the length of crack required to 

produce the same effect as an infinite strip. 
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The difficulties encountered in exj)erimentally holding a constant 
stress at the boundary as the crack moves is great enough to piirsuade 
most investigators to use displacement boundary conditions . 



Displacement Boundary Conditions 

For the same strip geometry subjected to displacement boiindary 
conditions on 

y = ±b u = 0 

Guess a displacement field that is everywhere zero for u and is odd 
in y for v. 

u(x,y) = 0 

O 

v(y) = + a^^y + a^y^ ... 

On the centerline of the sheet, v must be zero by symmetry; therefore 

a^ = 0. The strain field corresponding to these displacements is 

2 

e = 0 e = a^ + 3a_y + . . . e = e (x.y) Y = 

X y 1 3 z ' xy 



The stresses that resiilt from such boundary conditions are calciolated 
from the stress-strain law, 

0=0 - vd 

X y 



2 1 

a^ Sa^y + . . . = " vdjj] 



e_ = ^o„ + o 1 



z E x y 



From (B.lla) 



O = VO , 
X y’ 



and 



=1 + 3ajy2 ... =[^} 



(B.8) 

(B.9a) 

(B.9b) 

I 

(B.IO) 

(B.lla) 

(B.llb) 

(B.llc) 

(B.12a) 

(B.12b) 
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(B.12c) 



e = + v)<7 

z E' V 



The two dimensional equilibrium equations sin 5 >lify to 



da 



da 

r^ = o 



Using equation (B.12b), eqxiilibrium is satisfied if 



E 



5 [ 6 a y + ...] =0 

1 -v 

Setting a^ and all coefficients of higher order terms to zero, the 
solution reduces to 

u(x,y) = 0 

'^(y) = a^y 

V 

The boundary conditions are satisfied if a^ = therefore. 



u(x,y) = 0 



v(y) =^y 



The strains are obtained by differentiation. 

V 



e =0 

X 



e = ^ e = -v(l-v)r^ Y =0 
y b z ' 'b 'xy 



(B.13) 



(B.lUa) 

(B.l4b) 

(B.15) 



(b.i 6 ) 



and the stress field is expressed by 





(B.17) 



As before, when a crack is inserted into this specimen the constant 
stress field of equation (B. 17 ) is present at a distance out in front 
of the crack. As the crack advances, the domain of constant stress 
far in front of the crack is reduced by Ac, and the domain of zero 
stress far to the rear of the crack is increased by A 3 . The energy lost 
in the process of the crack moving by Ac reduces to the same expression 
as before, even though the stress and strain fields are different. 
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(B.18) 



AU = ^ a e 2bt Ac 
2 y y 

In the limit of infinitesimal crack increments 



dC 





2 



bt 



Once again from Irwins 



relationship in equation (B.6) 




or from (B.l?) 

K = a Vl-v^ 

1 y ^ 



(B.19) 



(B.20a) 



(B.20b) 
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